Gauge singlet scalar as inflaton and thermal relic dark matter 
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We show that, by adding a gauge singlet scalar 5 to the Standard Model which is non-minimally coupled 
to gravity, 5 can act both as the inflaton and as thermal relic dark matter. We obtain the allowed region of 
the (nis, ni/,) parameter space which gives a spectral index in agreement with observational bounds and also 
produces the observed dark matter density while not violating vacuum stability or non-perturbativity constraints. 
We show that, in contrast to the case of Higgs inflation, once quantum corrections are included the spectral index 
is significantly larger than the classical value (n = 0.966 for A' = 60) for all allowed values of the Higgs mass 
m/,. The range of Higgs mass compatible with the constraints is 145 GeV ^ni/, < 170 GeV. The 5 mass lies in 
the range 45 GeV < m,, < 1 TeV for the case of a real S scalar with large quartic self-coupling A,.,, with a smaller 
upper bound for smaller A.j. A region of the parameter space is accessible to direct searches at the LHC via 
h SS, while future direct dark matter searches should be able to significantly constrain the model. 

PACS numbers: 98.80.Cq, 95.35.+d, 12.60.-1 



I. INTRODUCTION 

Although the Standard Model (SM) of particle physics re- 
produces experimental results well, it does not provide a 
mechanism for inflation, dark matter, baryogenesis or neutrino 
masses. There are many Beyond the Standard Model (BSM) 
theories which attempt to do this. Such BSM theories usually 
involve a new scale between the SM scale and the Planck scale 
and many additional particles, with the SM then viewed as the 
low energy remnant of some more complete theory, which is 
only valid up to some cut off A < Mp. 

An alternative philosophy is to add to the SM the minimal 
number of new fields that are needed to address these issues. 
One example is the vMSM which is the SM extended 
by three singlet fermions to account for neutrino masses. In 
this case dark matter can be explained by a keV-scale ster- 
ile neutrino while baryogenesis occurs via leptogenesis due to 
sterile neutrino oscillations 101 • Therefore neutrino masses, 
dark matter and baryogenesis can all be explained within a 
very minimal extension of the SM, although this imposes non- 
trivial conditions on the sterile neutrino masses and couplings 
dlt]. A scale-invariant but very weakly-coupled scalar may 
also be added to serve as the inflaton |4]. Other minimal ex- 
tensions of the SM include the 'new minimal Standard Model' 
[S] and the 'minimal non-minimal Standard Model' |6]. One 
motivation for considering weak-scale extensions of the SM 
is the idea that the hierarchy problems of non-supersymmetric 
particle theories can be avoided if there is only one mass scale 
in the effective field theory below the Planck scale j^. 

Recently it has been suggested that inflation might be ex- 
plained purely within the framework of the SM, with the 
Higgs field itself serving as the inflaton (jstl . This has been ex- 
tensively investigated in a number of papers 
[l3, 14, 15, 16]. This is possible if the Higgs has a large non- 
minimal coupling to gravity. However, in order to account for 
dark matter, baryogenesis and neutrino masses, it is still nec- 
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essary to extend the SM. This might be achieved by combining 
Higgs inflation with the vMSM, but other extensions which 
are consistent with entirely weak scale particle physics could 
also be considered. In particular, it is well-known that stable 
particles with weak scale masses and electroweak strength in- 
teractions (WIMPs) produce a thermal relic density of dark 
matter which is naturally of the correct order of magnitude. 
Therefore there is a strong motivation to extend the SM by the 
addition of a particle with these properties. 

The aim of this paper is to propose an alternative 
minimally-extended version of the SM which is able to ex- 
plain both the mechanism for inflation and the presence of 
thermal relic dark matter. To this end, we add a stable gauge 
singlet scalar S to the SM. This is the simplest extension 
which obe ys gauge symmetry and can account for dark matter 
E [Sill IS Eh- A discrete Z2 or a global symmetry t/(l) 
must be imposed to ensure stability of the scalars; in the for- 
mer case it is natural to consider real scalars, in the latter case 
complex scalars. We then consider whether S can serve si- 
multaneously as a thermal relic dark matter particle and as the 
inflaton, producing the correct density of dark matter while 
at the same time obeying the observational constraints on the 
spectral index n and other inflation observables. Effectively 
we are replacing the Higgs scalar of Higgs inflation by the 
dark matter scalar S. As we will show, the model has the 
potential to relate particle physics, dark matter detection ex- 
periments and inflation observables, a connection that will be 
brought into focus in the near future by the LHC, the Planck 
satellite and future dark matter detectors. 

During the development of this paper a closely related 
model was proposed in ll22ll . This considers the same gauge 
singlet scalar extension of the SM to account for dark matter, 
but focuses on the case of Higgs inflation. As we will dis- 
cuss, there are some differences in the results for pure Higgs 
inflation, the model of |22] and our model, such that it may 
be within the reach of imminent experiments (Planck, LHC) 
to rule out or distinguish between Higgs inflation models and 
S'-inflation. 

Our paper is organized as follows. In Sec. |II]we introduce 
our model. In Sec. |III]we review the approach to calculat- 
ing radiative corrections in this class of models and derive the 
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renormalisation group (RG) equations. In Sec. |IV] we dis- 
cuss constraints coming from stability and perturbativity of 
the potential and slow-roll inflation observables. In Sec. IVlwe 
discuss S as dark matter, relating Xi,s and m,. In Sec. I VII we 
present our results and in Sec. I VIII we discuss our conclusions. 
Details of the derivation of the RG equations and the calcula- 
tion of the dark matter density are given in the appendices. 



II. THE 5-INFLATION MODEL 
A. Jordan and Einstein Frames 

The Jordan frame is the 'real world' frame, where we make 
measurements in a standard manner. The Einstein frame is re- 
lated to this by a conformal transformation which transforms 
the metric (and hence all other quantities) in a field dependent 
way. The usefulness of transforming to the Einstein frame is 
that it transforms away the non-minimal coupling to gravity, 
leaving the Lagrangian in a familiar form, where methods for 
calculating physical quantities are well known. What we think 
of as a conformal transformation to the Einstein frame is ac- 
tually composed of two separate parts. The first is a change 
of conformal frame, warping the metric, and the second re- 
defines the fields in a convenient form. Useful discussions of 
conformal transformations are given in 1I23I1 and lf24ll . 

Our procedure is to define the theory, including all radiative 
corrections, in the Jordan frame. We then transform to the Ein- 
stein frame in order to calculate the spectral index n, tensor- 
to-scalar ratio r and running of the spectral index a. The two 
frames are equivalent at low values of the fields. Since the in- 
flation observables are calculated when perturbations re-enter 
the horizon i.e. at late times when the fields are small, the 
results calculated in the Einstein frame are the same as if we 
had calculated them with the non-minimally coupled scalar 
field in the Jordan frame. 



^ is the constant contribution to the total S mass squared, 
m^, which also gains a contribution from the coupling to the 
Higgs. For now we consider only the physical Higgs field h, 

V2 [ h + 



where H 



and h is real. We choose the direc- 



tion of inflation such that S — where s is real. 

Our aim is to calculate the inflation observables n, r and a. 
This is best done using the slow-roll approximation, which 
cannot easily be formulated in the Jordan frame. We will 
therefore make a transformation of the whole action, includ- 
ing radiative corrections, to the Einstein frame, redefining the 
fields (s — > Xs, h — > X/7) to ensure canonical normalisation. We 
then compute the slow roll parameters in the Einstein frame, 
using the coupling constants which we have run (in the Jor- 
dan frame) to the appropriate scale. Quantities in the Einstein 
frame will be denoted by a tilde (e.g. ^^v)- From here on we 
set M = Mp (reduced Planck mass), since the correction to M 
due to the Higgs expectation value is tiny compared with Mp. 

For general h and s the transformation to the Einstein frame 
is defined by 



with 



£1^ 



Ml 



Ml 



(3) 



(4) 



The fields are redefined by 

ds 

and 



la^+e^y/Ml 

£24 



(5) 



dh y 



(6) 



B. Non-minimally Coupled Gauge Singlet Scalar Extension of 
tlieSM 

We define the action in the Jordan frame to be 



Sj = J V^d'x[^LsM+{d,Hnd^-H) + {d,Snd^S) 

^ ^irH'HR - ^,5^5/; - ViS^S,H^H)j (1) 

where V{S'^S,H^H) = y W + y(i) + • • .. Here 



resulting in the Einstein frame action 



mIr U 



where 



-SM 2 2 



MlQ.* dx»dXh 



(7) 



(8) 



„2\ ^ 



(2) 



is the tree-level potential and V^^\V^^\ ... are the 1-loop and 
higher-order quantum corrections. L-^ is the Standard Model 
Lagrangian density minus the purely Higgs doublet terms. 



and 



U{Xs,Xh) = ^vi''^) 



with 

{9) 
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We will be interested in inflation purely along the s direc- 
tion. (Inflation in the h direction for real S was considered in 

m.) In this case, h = 0, A{x,,Xh) = and £2^ ^ 1 + 

For s ^ Mpj \f\, which is relevant for inflation, the classical 
potential in the Einstein frame becomes 



-exp 



-2 



(10) 



This is shown in Fig. [T] Thus U{%s,Q) °^ 1/^?- Similarly, 



U(Zs, 0) / M, 




FIG. 1 : Classical potential in Einstein frame, in limit s S> Mp / 'Ji^. This 
figure is plotted for real S with hi/, = 160 GeV and Xs = 0-2. Inflation occurs 
along the exponentially flat plateau. 

along the h direction with 5 = 0, U{0,Xi,) °^ 1/^^- Therefore 
if t,s ^ t,h, the minimum of the potential at large s and h will 
be very close the h = direction and so inflation will naturally 
occur along the s direction. In the following we will consider 
the limit where the Higgs boson is minimally coupled to the 
Ricci scalar at the weak scale, = 0, but we allow for its 
running by including the RG equation for . 



III. RADIATIVE CORRECTIONS 

Our strategy is to calculate quantum corrections to the tree 
level potential in the Jordan frame. To do this we use the RG 
equations to run the couplings from the SM scale to the infla- 
tion scale. We then use those values of the coupling constants 
to calculate the Coleman- Weinberg correction to the potential, 
V ^ yf") + miM], where V is given by Eq. ©. This 
is then transformed to the Einstein frame to study slow-roll 
inflation. 



A. Coleman-Weinberg Potential 

Constraints on the scalar couplings will come from the sta- 
bility of the electroweak vacuum and the requirement that the 
potential remains perturbative for field values less than Mp. 
We therefore impose the conditions for vacuum stability and 
perturbativity along both the h = Q and 5 = directions. To do 



this we derive the Coleman-Weinberg potential for each direc- 
tion (s or h). The i-direction potential is used in the slow-roll 
inflation calculations, while both i-direction and /i-direction 
potentials are used to calculate the stability of the potential 
and to check perturbativity. We use the MS renormalisation 
scheme throughout. The one-loop potential for the s direction 
in the MS scheme is ll26ll 



-H? In- 



4^-* \ 2 
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where 



^ [m'^^+XsS (complex 5). 



The one-loop correction for the h direction is 
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where 
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Ph 
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ml + 3ci,Xhh^, G, 



2 , ,2 ^,^„,2^X,/22, 



1 



CsXhsh and 



(real S) 

- jkhsh^ (complex S). 



(14) 



In these equations q and c/, are suppression factors to be dis- 
cussed below. 



B. Suppression of Scalar Propagators 

As the fields appearing in the RG equations are quantised 
in the Jordan frame, the commutation relation for an arbitrary 
scalar (|) 



is satisfied, where ^£ 



n = 



dL 




(15) 



(16) 
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and r|^ = (1,0,0,0). This is obtained by transforming the ac- 
tion to the Einstein frame but not redefining the scalar fields 
lfl4ll . Inserting Eq. (fTST i into Eq. (flST l and writing in terms of 
the Jordan frame metric gives 



and so. 



where c((|)) 



a2 



Therefore the commutator and hence 



the scalar propagator will be suppressed by a factor c 

^ 1 



In the case of minimally coupled scalars, ^3^^ 

c((|)) = 1. In our case, both s and h are in principle sup- 
pressed by c{s) and c{h) respectively. In practice, we set either 
c/,(= c{h)) = 1 or Cj(= c{s)) = 1, depending on the direction 
of the potential being considered. The suppression factor C(|, 
(where (|) is i or h) is then 



C(|) 



1 + (6^.^ + 1) 



(18) 



When calculating the RG equations or Coleman-Weinberg po- 
tential, one suppression factor is inserted for each /; or s prop- 
agator in a loop but not for the scalars corresponding to imag- 
inary part of S or the unphysical degrees of freedom of H. 
The suppression factors will have a significant effect on the 
running of the scalar couplings. 

In the context of Higgs inflation, it was shown in iEtI Izsll 
that unitarity breaks down in tree-level Higgs-graviton scat- 
tering processes at energies E ~ Mp/t,/,, due to the large non- 
minimal coupling to gravity. There are two possible ways to 
interpret this. In [Jj]] and u3l it is suggested that the appar- 
ent breakdown of unitarity should be interpreted as a change 
in the nature of the Higgs degree of freedom, rather than as a 
cut-off for new physics. In [13] it was observed that at field 
strengths (h^) ^ Mp /^|, the Higgs scalar h no longer behaves 
as a canonically normalized scalar, resulting in suppression 
of the Higgs propagator as discussed above. The results for 
Higgs scattering cross-sections at the corresponding energies 
are therefore expected to be modified, which may justify the 
extension of the theory to energies greater than Mp/t,i,. In 
ifTTll . it is proposed that the onset of unitarity violation could 
indicate a change in the dynamics of the Standard Model to a 
strongly coupled regime. This is described by the 'chiral elec- 
troweak theory', which is equivalent to the Standard Model 
with the radial Higgs degree of freedom frozen'. 



' Since the chiral electroweak theory has the radial Higgs mode frozen, it 
explicitly breaks SU (2)l x 6'( 1 )y symmetry, making it non-renormalizable 
IIILI22I . We note that explicit gauge symmetry breaking is avoided when 
a gauge singlet scalar plays the role of the inflaton. 



A more conservative point of view is to restrict the model 
to the regime where the semi-classical and adiabatic ap- 
proximations are both valid iEtIi . These approximations are 
vahd when 1 > ///M > where H w VXJ^Mp/^i and 
M "Z Mp/i^s is the scale where new physics becomes impor- 
tant. As A-s in our model is not constrained by phenomenol- 
ogy (unlike the case of Higgs inflation where A,/, is fixed by 
the Higgs mass), it would be possible to have <C 1 . 



C. Initial Conditions 

We take the initial values of the coupling constants to be 
defined at the renormalisation scale /j = m,, with m, = 171.0 
GeV and v = 246.22 GeV. The gauge couplings are given by 



g {mt) 

gl{mt) 
4n 



= 0.03344, ^-4^ = 0.01027 and 



4k 



= 0.1071. 



(19) 



g and g' are obtained by an RG flow from their values at 
/J — Mz, which are given in ll29ll . while gs is calculated nu- 
merically. (See 1 1 3l1 and references within for details.) 

We use the pole mass matching scheme for A,/,(m,) and 
yt{mt) as detailed in the appendix of lf30ll . This relates the 
physical pole masses to the couplings in the MS renormal- 
ization scheme. The remaining coupling constants are not 
fixed by observation and we are free to choose them. We take 
^h{fnt) = and choose ^.v('«r) such that the model is correctly 
normahsed to the COBE results at the inflation scale li3li,i3Z1 : 



- = (0.0027 IMp)^ 



A,s(m,) is not directly measurable and so we take two reason- 
able values: 0.2 and 0.025. The higher of these corresponds to 
Xs{mt) close to the perturbativity limit. A-/,s(m,) is treated as a 
free parameter although it is in principle measurable through 
the thermal relic S dark matter density and scattering rate in 
dark matter detectors, as well as through the Higgs decay 
width to S pairs should it be kinematically possible. 



D. Renormalisation Group Equations 

In our analysis we use the two-loop RG equations for the 
SM and modify these to include the leading order contribu- 
tions of S. We also include the propagator suppression factors 
for the s and h directions. We refer the reader to 1 14] and |3(3] 
for the SM one and two loop equations (including only the t- 
quark Yukawa coupling), reproducing here only those which 
are modified by the addition of the S particle. Using the tech- 
nique detailed in [33] and further discussed in Appendix A, 
we find for the one-loop p-functions of the scalar couplings. 



16^^Pl!.' 



5 



+ (- 



3/ + 12y2)X/, 
(real S) 



■hs 



I A,^^. (complex S), 



(20) 



) ?^/,.v + 6^, A./,, 

( 6c'^Xshrs (real 5) 

'\ (6c2 + 2)A.,A./,, (complex 5) ^ ^ 



and 



167t^P 



.2r(1) 



(4 + 3)^1 



18c2X2 

fl8c? + 



(real 5) 
I A,^ (complex S), 



(22) 



-7-, ? = In-^ and u is the renormalisation scale. 

at ' " 



where h^§, 

We choose the value of /j in order to keep the log terms in the 
Coleman- Weinberg potential small, setting /j — S(,o, where S(,o 
is the field value 60 e-foldings before the end of inflation. In 
Appendix A we relate the gauge singlet model to the matrices 
defined in [33] which are used to compute the RG equations. 

We also obtained the RG equations for the non-minimal 
couplings to one-loop. The details of this calculation are given 
in Appendix A. The resulting equations are 



1671- 



dt 



(3 + c/,)A/„ ( 



1 

6 

6C.sA; 



(real S) 



(6Ci + 2) "kg (complex S) 



,(23) 




FIG. 2: Running of scalar couplings showing the effect of suppressing the 
i propagator. The dash line shows the coiTesponding value of when /.i = 
s. This figure is plotted for real S with m;, = 160 GeV, ^^(m,) = 0.2 and 
'k,^(m,) = 0.\. 
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FIG. 3: Running of the non-minimal couplings of i and h to the Ricci scalar. 
In this, is set to zero at |U = m, . It can be seen that ^ throughout. The 
figure is plotted for real 5 with m;, = 160 GeV, Xs{m,) = 0.2 and Xi,^{mr) = 
0.1. 



and 



16r 



dt 



(6 + 6c/,) A/, 



-63'?--(3g'+g' 



1 

c.Xhs (real S) 

( 1 + Cj) A/,j (complex S) 



(24) 



We show in Fig. |2] the running of the scalar coupling con- 
stants A/,, A.V and A/,., in the s direction for the case of a real 
S dark matter particle, assuming a small value of Xs{mt) = 
0.025. The figures are plotted in terms of f = In/j/mr, from 
/J = nit fj — Mp. We also show the suppression factor Cj, 
Eq. ( fTSl ), to demonstrate its effect on the running of the cou- 
plings. In Fig. [3] we plot the running of t,^ and for the 
same initial conditions as in Fig. |2] We observe increas- 
ing from its initial value of at f = 0, but always remaining 
much smaller than t,^. This is important for the consistency of 
our model since inflation will occur along the s direction only 
if t,s ^ ^/j. Otherwise inflation would be expected to occur 
along a more general flat direction in the (s, h) plane. 



IV. CONSTRAINTS 

We calculate the bounds on m/, and Xhs{mt) by applying 
three constraints: (i) stability of the electroweak vacuum, (ii) 
perturbativity of the potential and (iii) consistency with the 
observed spectral index n and with limits on the tensor-to- 
scalar ratio r and running spectral index a. A possible fourth 
constraint, 'wrong-way-rolF (^ > 0), which plays a role in 
Higgs inflation L221 . is generally not violated in our model. 



A. Vacuum stability and perturbativity 

We require stability of the electroweak vacuum for s and h 
up to Mp. (We do not consider the possibility of a metastable 
vacuum, which depends on the cosmological evolution of the 
vacuum state.) This imposes the constraints A; > 0,A/, > 
and either A/jj > or A^^. < 4A/,Ai. We will check the stability 
of the vacuum in both the s and the h direction. In practice this 
means that we run the RG equations with c, = 1 and varying 



6 



c/, and again with c/, = 1 and varying c^. 

We also require the couphng constants to He within the per- 
turbative regime up to the Planck scale in both the s direc- 
tion and the h direction. We apply the perturbativity condition 
A,J < 4n to the coupling constants X'j defined through the po- 
tential 



1 



(25) 



The couplings in this potential appear in the Feynman vertices 
without additional numerical factors. < 471 then ensures 
that loop corrections are smaller than tree-level processes. 
(Above 471, the coupling constants quickly grow towards a 
Landau pole. Therefore altering the definition of perturba- 
tivity will not significantly change our results.) This leads 
to the conditions on the couplings as defined in our potential 
A-/,, Xs <2k/3 and A,/„ < 47i. 



where the end of inflation is defined by fj = 1 . Although this is 
calculated in the Einstein frame, it is straightforward to show 
that .^V is equal to the number of e-foldings in the Jordan frame 

up to a small coiTection"*, N w + ln(l/VA^). We will 
use .^V = 60 when calculating inflation observables. This is a 
reasonable assumption given that the reheating temperature in 
this model will be high^. 

Using the tree-level potential and the approximation 



1 we estimate the tree-level slow-roll parameters 



to be e 

f2 Al 



Mi 



andC^^^I^, where 4 



4 '"p ~ 4"'p 

4-MpN/3^s- A calculation of the tree-level spectral index, 
tensor-to-scalar ratio and running spectral index gives 



,(0), 



O 



1 



= 0.966 ; 



B. Constraints from slow-roll inflation 

The present observational constraints on inflation are n = 
0.960 ±0.013 (1-a), r < 0.22 and -0.068 < a < 0.012 jsl . 
Inflation occurs through the standard slow-roll mechanism, 
which we formulate in the Einstein frame. The potential in 
the direction is 

where U^^\s) is given by Eq. ( fTTT ). The slow roll parameters 

2 

are 

^_M^/i^y 



11 = 



2 \UdxsJ ' 
U dxl 



and 



_ M^dU_fO_ 



f/2 dXs dxl ' 

From these we can calculate the observable quantities 

n = 1 -6e + 2fi , 
r = 16e 



(26) 



and 



a = 



dn 
dink 



= -16iie + 2462 + 2^2 



(27) 



The number of e-foldings of inflation is given by the standard 
expression jsTIl 



1 U 



(28) 



Xend p 



dXS 



■ We use ^ rather than 4 to avoid confusion. 



r = 3.3x 10 



-3 , 



a = 6.2x 10 



-4 



Thus r and a are negligibly small when compared with the 
observational limits. 

Radiative corrections have a significant effect on the slow- 
roll parameters. This is not surprising, as the tree level po- 
tential is exponentially flat and the radiative corrections add a 
small but significant slope. Including radiative corrections we 
find 



Mi 



ds 
dXs 



and 



4 

-3 a 



FX 

s 



3,2 



Mj 



FX 



where 



X 
F 



(l+Flns + oy 



1 



1 



wx:^4^''^^ 

9c]X] 



3) XL 



(real S) 



and 



D 



{9c] + \)X] (complex 5) 



^IXl 



(29) 



(30) 



(31) 



(32) 



\6%^Xs 



9c2x2ln3£A 
9^212 In ^ 



^ln£^!^ + ^ln^- 
4 2^P- 4 2^P- 

(real S) \ 

A,2ln% (complex 5). J 



(33) 



^ yV is defined as ln(fl„,rf /cin), where a is the scale factor of the Friedmann- 
Robertson-Walker metric in the Einstein frame, which is related to the con- 
ventional scale factor by a = Q.a. 

^ We note that N = 55 gives n = 0.963, therefore n is not particularly sen- 
sitive to the value of N. We plan to compute the reheating temperature 
precisely in a future paper, which will fix N. 
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The terms originating from U^^^ are subdominant in fj, but 
for a range of values of A,/,j and A,, they can become more 
important than the tree-level result in 8. 

We calculate the field value at 60 e-foldings before the end 
of inflation as follows. First we calculate Send using jfjl = 1. 
This gives (tree level) 



VI. RESULTS 



4Ml 
3 ^. 



(34) 



Then the standard expression Eq. ( |28] ) is integrated using Eq. 
I and the approximation X = constant to give 



N ~ yln 



4+ran2(w)y 4"lv^i'(w) 



with £i} as defined in Eq. (|4]i and 



(35) 



(36) 



V. THERMAL RELIC DARK MATTER 

We assume that dark matter is due to thermal relic gauge 
singlet scalars. The non-minimal coupling to gravity will 
not affect the S dark matter density as the field is at very 
low values compared to Mp. If we assume that a gauge sin- 
glet scalar is responsible for the observed dark matter den- 
sity iO-Duh^ = Pdm/Pc- ^ 0.1131 ±0.0034 H) then we ob- 
tain a relationship between njj and Xhs{mt). We use the Lee- 
Weinberg approximation |35] to calculate the relic density of 
S. This is discussed in Appendix B, where we also review the 
annihilation cross-sections. For a given A-/,j and m/, there are 
up to four corresponding values of m^ . An example is shown 
in Fig.|4]for the case where m/, = 160 GeV. The cusp-like fea- 
ture is due to S annihilations to WW and ZZ pairs close to the 
Higgs pole. In this region the S mass is relatively insensitive 
to Xhs- Note also that large values of A,/,., are possible for 
slightly below the Higgs pole. 
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FIG. 4: The value of as a function of X;,,(mf) necessary to produce the 
correct density of thermal relic dark matter. In this example m;, = 160.0 GeV. 
The solid line indicates real S and the dashed line complex S scalars. 



InFig. [5(a)] we show the case of real S with 'small' Xs{mt) — 
0.025. The range of allowed Higgs mass is 145 GeV ^rrih^ 
170 GeV, where the lower bound is from vacuum stability 
in the h direction combined with 5-year WMAP l-o upper 
bound n < 0.973 and the upper bound is from perturbativity 
of Xf, in the s direction. The corresponding range of Xhs{mt) 
is \khs{fnt)\ ^ 0.15. Larger values of n allow larger \khs{mt)\, 
up to an upper bound |A,/,((m,)| sa 0.55 (at n > 0.980), which 
comes from the perturbativity bound on in the h direction. 
In this case the lower bound on the allowed Higgs masses is 
shifted downwards to 130 GeV < m/, < 170 GeV. In Fig. [5(b)l 
we show the corresponding results for complex S. The al- 
lowed parameter space is very similar to the case of real S. 

In Fig. |5(c)| we show the results for the case of 'large' 
Xs{mt) — 0.2. In this case the range of Higgs mass is similar to 
the small Xs{mt) case, but now the origin of the bound is per- 
turbativity of Xs in the h direction rather than the WMAP up- 
per bound on n. As Xs{mt) increases from 0.2, the allowed pa- 
rameter space will rapidly diminish due to the decrease of the 
A,j perturbativity upper bound on X/„(m;). As seen in Fig. |5(d)[ 
the allowed parameter space vanishes for the corresponding 
case with complex S. 

As discussed in Sec. lIIIBl when Xs ^ 1 we can avoid po- 
tential problems due to unitarity violation. Choosing a small 
value of Xs at the scale of inflation to satisfy this bound should 
not be a problem. Comparing Fig. |5(a)| and Fig. |5(c)| we see 
the bound on Xhs due to n decreases with Xs, therefore we 
would expect the allowed range to decrease further with even 
smaller Xs- A smaller range of Xhs will tend to drive nis closer 
to the Higgs pole (Fig.|4|, increasing the chances that it could 
be detected in the near future. 

An important point is that the value of n can be significantly 
larger than the classical value n = 0.966 over the whole range 
of allowed Higgs mass. In Fig.|6]we show an example of the 
variation of n with m/, for fixed Xks and A,,. This contrasts with 
the case of Higgs inflation without additional scalars, where a 
significant increase of n relative to the classical value is pos- 
sible only for a small range of Higgs mass close to the vac- 
uum stability lower bound; in fl^] a significant increase of n 
from the classical value is obtained only for m/j 132 GeV. 
Therefore if Planck, which will measure n to a 2-o accuracy of 
±0.005, should find « significandy larger than 0.966 + 0.005 
while LHC finds a Higgs with mass larger than 135 GeV, then 
S inflation will be compatible with the observations but Higgs 
inflation will be ruled out. We may also compare the range of 
Higgs mass allowed by S inflation with that allowed by vac- 
uum stability and perturbativity in the Standard Model. In f3^ 
the range is given as 128.6 GeV < mi, < 175 GeV, where the 
lower bound is from vacuum stability and the upper bound is 
from perturbativity of the Higgs self-coupling up to Mp. We 
see that the allowed range in S inflation is somewhat narrower 
Therefore S inflation may be ruled out relative to the conven- 
tional SM if m/, is observed close to the SM lower or upper 
bound. 

Fig. 7 shows the range of iris and m/, which is consis- 
tent with S inflation and thermal relic S dark matter when 
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(a) Real 5, A,(m,) = 0.025 (b) Complex S, him,) = 0.025 




110 120 130 140 150 160 170 ^20 130 140 150 160 

(c) Real 5, X,(m,) = 0.2 (d) Complex 5, \,(m,) = 0.2 



FIG. 5: Allowed region for inflation in the .s-direction. Excluded regions are shown in grey. Limits from couplings in the .v-direction are shown with dashed 
lines, those from the couplings running in the /i-direction have solid lines and the 1-a upper limit on n is dot-dashed. In (a) we show the line n = 0.981 
(dot-dot-dash) demonstrating the variation of n. 



0.98 


WMAP + 1 t7 


X„,(m,) = 0.16/ 


0.97 


5 inflation classical value 




0.96 


WMAP central value 





110 120 130 140 150 160 170 



FIG. 6: Showing the variation of n with m;, for X,{mi) = 0.025 and ^;,,(mf) = 
0.16. WMAP central value and 1 a upper bound shown with short dashed 
lines; classical n for S inflation shown with dashed line. 



n < 0.973 and all vacuum stability and perturb ativity con- 
straints are satisfied. We also show the line = m/,/2, which 
is the limit at which it is possible to pair produce S scalars 
via Higgs decay at the LHC |37]. For the case of real S and 
'small' Xs{mt) — 0.025, Fig. |7(a)[ we see that nis is mostly in 



the range 50 GeV < < 500 GeV, reaching 750 GeV close 
to its lower bound. For complex S, is more constrained, 
with values in the range 50-500 GeV. This can be easily un- 
derstood since the dark matter density for a complex S is twice 
that for a real S of the same mass, therefore a smaller mass is 
required to produce the same density. From Fig. |7(c)| we see 
that a larger value of Xs{mt), Xs{nit) — 0.2, permits a wider 
range of S mass, with in the range 45 GeV to 1 TeV If 
we instead considered the 2-o WMAP bound, the parameter 
space in Fig. |7(a)| and Fig. |7(b)| would increase, while Fig. |7(c)| 
would be unchanged. 

We note that while a large region of the allowed parameter 
space is at values of the S mass which are large compared with 
the weak scale, there is no reason to expect the S mass to be so 
large. The S mass squared is — m^^ + A,/,jV^/2. Therefore if 
ms„ is of the order of the weak scale, which is the most natural 
possibility in a theory based on a single mass scale, we would 
expect nis to be no larger than a few hundred GeV. 

There is a small region of the parameter space which sat- 
isfies nis < OT/,/2, with the lower bound on in the allowed 
region being slightly below m/,/2. This means that it is pos- 
sible for the S inflaton to be produced at the LHC via Higgs 
decay [37]. Thermal relic S dark matter would then originate 
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800 



600 



400 



200 




1000 



800 



600 



400 



200 



(c) Real5, X,(m,) =0.2 
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(d) Complex 5, (m, ) =0.2 



FIG. 7: Allowed region for inflation in the ^-direction, with a l-o upper limit on n. Excluded regions are shown in grey and all masses are in GeV. The dashed 
line shows mi, = 2m,. Below this line, production of 5-particles at the LHC (via SS' decay) is possible. There is no allowed region in (d). 



from freeze-out of near resonant S annihilation to WW and ZZ 
close to the Higgs pole, slightly below the Higgs pole also 
implies that Xhs{mt) can be large, as can be seen from Fig. |4l 
Therefore if iris < 'w/,/2 then the 5-nucleon scattering cross- 
section due to Higgs exchange is likely to be large, enhancing 
the possibility of observing S dark matter in direct detection 
experiments. 

Collider and direct dark matter detection experiments 
should be able to constrain the allowed parameter space^. 
Combined data from the DO and CDF collaborations show that 
a Higgs boson mass in the range 160 GeV < m/, < 170 GeV is 
excluded at 95% confidence level |39]. This exclusion reduces 
the available parameter space of the model by a little less than 
a half, however it does not make a large difference to the range 
of nis, as this is largest at low values of m/,. Present bounds on 
direct detection of S dark matter from XENON 10 and CDM- 
SII rule out S mass in the range 10 GeV to (50,70,75) GeV for 
Higgs masses (120,200,350) GeV |'37','4^. Comparing with 
Fig. 7, we see that the upper bound from direct detection is 
already close to the lower bound on the range of nis allowed 



^ y-ray and antimatter signals can also constrain the model IbsIi . 



by the S'-inflation model. Thus although most of the param- 
eter space is allowed at present, a substantial part of the (nis, 
nil,) parameter space will be accessible to future dark matter 
detectors. 



VII. CONCLUSIONS 

We have shown that a gauge singlet scalar can serve simul- 
taneously as the inflaton and as a thermal relic dark matter 
particle. As in the case of Higgs inflation, this requires a 
specific large non-minimal coupling of S to gravity. Consis- 
tency of the model with (i) stability of the electroweak vac- 
uum, (ii) perturbativity of the scalar potential as a function 
of s and h up to the Planck scale and (iii) the observed spec- 
tral index, constrains the (m, ) , m/, ) parameter space. (The 
tensor-to-scalar ratio r and the running of the spectral index 
a are both negligibly small compared with the observational 
limits.) Imposing the 5-year WMAP l-o bound n < 0.973 
implies that the range of the coupling of S to the Higgs is 
|A./jj(m()| < 0.15, which can increase up to |A,/„.(m,)| < 0.55 
for small S self-coupling and larger n. The range of Higgs 
masses is similar to but not identical to that of the Stan- 
dard Model, with 145 GeV < m/, < 170 GeV for « < 0.973 
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and small X,{m,), shifting to 130 GeV < m/, < 170 GeV for 
n ^ 0.980. Demanding that the S annihilation rate through 
'Xhs{mt) produces the correct thermal relic S dark matter den- 
sity translates each Xi,s{mt) into a discrete set of possible val- 
ues of nis- Combined with the above constraints this deter- 
mines a range of nis and m/, which is simultaneously consis- 
tent with thermal relic dark matter and a stable, perturbative 
scalar potential which can account for the observed spectral 
index. The range of is sensitive to and to whether 

S is real or complex; for Xs{m,) — 0.025 and real S the range 
is 50 GeV ^ ^ 750 GeV, with the upper limit increasing 
to 1 TeV for Xs{mt) = 0.2. For complex S the range of is 
naiTower, 50 GeV < m, < 500 GeV for A.,(m,) = 0.025. 

Comparing with Higgs inflation in the unextended SM, a 
key difference is the range of possible n versus m/,. In the 
Higgs inflation case this is expected to be very close to the 
classical value n — 0.966 except for mi, close to the vacuum 
stability limit, m/, 130 GeV. In S inflation a significant de- 
viation from the classical value is expected over the whole 
range of m/,. Planck is expected to be able to observe n to 
an accuracy of ±0.005 (2-0). Therefore Planck could provide 
evidence in favour of S inflation relative to Higgs inflation, 
depending on what value of mh is observed at the LHC. 

Comparing with the model of f23\, which is based on the 
same gauge singlet dark matter model but considers inflation 
along the Higgs direction, a notable difference is that in our 
model n is strictly larger than the classical value n — 0.966, 
whereas in ||22|] it is possible for the spectral index to be- 
come smaller than the classical value. In addition, the results 
of 1 22] indicate that the spectral index becomes close to or 
smaller than the classical value at rrif, ^ 160 GeV (see Fig. 6 
and 7 of (tS]), whereas in our model the deviation from the 
classical value becomes larger as m/, increases. Therefore it 
may be possible to distinguish between 5-inflation and infla- 
tion along the Higgs direction, depending on the Higgs mass 
and the spectral index. 

There is a small region of the (m, , m/,) parameter space with 
nis close to m/,/2 which is consistent with production of S par- 
ticles at the LHC via Higgs decay, h SS. Therefore if S is 
observed at the LHC then thermal S dark matter must origi- 
nate from freeze-out of near resonant S annihilations to WW 
and ZZ at the Higgs pole. This also allows Xi,s{m,) to be large 
and so a significantly large Higgs decay branching ratio may 
be expected. A large 5-nucleon scattering cross-section via 
Higgs exchange is also expected in this case, which should 
allow the parameter space of the model to be probed by fu- 
ture dark matter detection experiments. The lower bound on 
nis is close to the present upper bound on nis from direct dark 
matter detectors (XENONIO, CDMSll). Therefore significant 
constraints on the parameter space (or possibly direct detec- 
tion of S dark matter) may be expected in the future as dark 
matter detectors improve in sensitivity. 

It is natural to ask how predictive the 5-inflation model can 
be. The barrier to a precisely predictive model is the depen- 
dence on the S self-coupling Xs{m,), which is not directly ob- 
servable. In principle, there are four observable quantities: «, 
nil,, nis and Xks{nit)- The input parameters of the S'-inflation 
model are m/,, nis, Xhs{nit) and Xs{mt). Therefore n cannot 



be predicted exactly as there will always be a dependence on 
Xs{nit), even if the other parameters of the model are fixed 
by experiment. Nevertheless, as we have shown, the possible 
range of n can be constrained by vacuum stability and pertur- 
bativity constraints which constrain Xs{nir). In addition, in the 
limit of small Xs {nit ) the model could become effectively in- 
dependent of Xs{ni,). In this case we can predict n if ;«/,, nis 
and Xhs{nit) are fixed by the LHC and direct dark matter ex- 
periments. If we are fortunate enough that S'-inflation occurs 
in this limit, then the model can be predictive and testable. 

We have focused on the case of inflation along the s direc- 
tion. In general, inflation could occur along a more general flat 
direction in the is, h) plane, depending on the value of 
and the V{s,h) scalar potential couplings. 5-inflation may be 
expected to be a good approximation in the limit ^ ^/,, as- 
suming that the couplings Xs, Xf, and A,/,, are of a similar order 
of magnitude. The opposite limit of inflation along the h direc- 
tion was studied in ll22ll . This may be expected if ^ ^s- The 
case of inflation along a general trajectory in the s, h plane, 
which would be expected if and are the same order of 
magnitude, remains to be investigated. 

The 5-inflation model provides a model for inflation and 
dark matter which is based purely on weak scale particles and 
interactions. In order to have a complete model of cosmology, 
we also need to address the issues of reheating and the origin 
of the baryon asymmetry. Reheating will be very similar to 
the case of Higgs inflation. In Higgs inflation reheating oc- 
curs via parametric resonance of the oscillating Higgs field to 
W bosons via the interaction Lgj, J.5]. In a similar 

way, in 5-inflation reheating will occur via parametric reso- 
nance of S oscillations to Higgs bosons via the |5|2|//|2 inter- 
action. Baryogenesis could occur via the oscillating leptoge- 
nesis mechanism 101 or low-scale resonant leptogenesis ll4Tll 
once the SM is extended by sterile neutrinos in order to ac- 
count for neutrino masses. Alternatively, baryogenesis could 
occur via electroweak baryogenesis, which may be possible 
in scalar extensions of the SM. Additional scalars interacting 
with the Higgs can produce a sufficiently strong 1st order elec- 
troweak phase transition. This usually requires that the gauge 
singlet scalar gains a vacuum expectation value (vev) after the 
transition ll42ll . therefore a more complicated model with two 
or more additional scalars would be required ^. 

In contrast to many inflation models, S'-inflation is notable 
for the close relationship it implies between the observables 
of inflation (in particular, the spectral index), particle physics 
(in particular, the Higgs mass and Higgs decay width) and the 
direct detection of dark matter. It can therefore be directly 
tested by the experimental and observational advances which 
are anticipated in the near future as the LHC, Planck satellite 
and future direct dark matter detection experiments come to 
fruition. 



* We note that it may be possible to evade this if the scalar has an expecta- 
tion value prior to and during the electroweak phase transition but its vev 
vanishes in the vacuum after the transition H8i1 



11 



Acknowledgements 



The U{1)y generator is 



This work was supported by the European Union through 
the Marie Curie Research and Training Network "Uni- 
verseNet" (MRTN-CT-2006-035863). The work of R.L. was 
also supported by STFC. 



Appendix A: Renormalization Group Equations 

In this appendix we relate the gauge singlet model to the 
notation of [33] used to compute the RG equations. We also 
briefly review the computation of the RG equations for t,^ and 



A. RG equations for scalar couplings 

In ll33ll the general RG equations are given to two-loops in 
the MS scheme. The anomalous dimensions and p-functions 
are expressed in terms of real (reducible) representations of 
the scalar fields and Majorana spinors. To obtain the modifi- 
cation of the RG equations due to the S field, we express the 
Higgs doublet and gauge singlet scalars as a set of six real 
scalar fields, (|), (/ = 1...6), where 



H 



J_ / (j)l +/(j)2 



and 



(A-1) 



(A-2) 



(For the case of real S, (|)6 = 0.) 

Writing the Higgs doublet as a real representation in the 
form ((|)i,(|)2,(|)3,(|)4)^, the SU{2)l generators (9^^ in the nota- 
tion of 1330) are 
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0/00 
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(A-4) 



and 
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0/0 



(A-5) 



F 
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F 

0-70 



(A-6) 



where F = 1/2 is the hypercharge of the complex fields in the 
Higgs doublet. 

The only Yukawa coupling we consider is the top quark 
Yukawa coupling. In 4-component spinor notation this is (in 
the notation of 1331] ) 



qH'^^'^'q + h.c. , 



(A-7) 



where H_ is the Yukawa coupling matrix, q ~ {ui,di)^ is the 
SU{2)i quark doublet and (|) is the Higgs doublet. In our case 



^//(j)''' "q = tRy,tL(if° - tRy,bL(i;>^ 



(A-8) 



(In this we have suppressed colour indices.) We define 
a reducible representation (in the notation of (y^) by 
(\(/i,\(/2,V(/3) = {t^,tL,bL), where ti, bi and ai-e the two- 
component spinors which form the Dirac spinors in the chi- 
ral representation (t = {ti^tRY etc), with = — /02f^. The 
Yukawa coupling can then be written as 



y/j\(/,^\(/^(j)„ + /!.c. (fl= 1,2,3,4) 



where 





(A-9) 



(A-10) 



(A-11) 



(A-12) 



and 




(A-13) 



The corresponding SU{2)l generators acting on i|/ are 





fi = i I 1 
^ * 1 




(A- 14) 



(A-15) 
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and 




(A-16) 



The U{1)y generator is 




(A-17) 



(Suppressed colour indices should be summed over when tak- 
ing traces in the formulae of ifssll .) Finally, K = 1/2 should 
be used in fi3\ since are two-component spinors. With 
these definitions of 9^, Y" and t^, the formulae in [fH can 
be used to compute the RG equations to two-loop order as a 
function of the t-quark Yukawa coupling, gauge couplings and 
the scalar couplings. 



density. We will approximate (av,.^/) by the centre-of-mass 
cross-section for non-relativistic S annihilation. The cross- 
sections for real and complex S are the same; we will present 
results for the real case. The tree-level processes contributing 
to S annihilation are (i) SS hh, (ii) SS WW, (iii) SS ZZ 
and (iv) SS ff (where / is a Standard Model fermion). (i) 
proceeds via a 4-point contact interaction, an s-channel Higgs 
exchange interaction and a t- and u-channel S exchange inter- 
action. The resulting (ovi-ei) is 



{aVrel)hh = 



647tOT? 



3ml 



2X 



■hsV 



2 

mf 
mi 



{Am^ — m 

1/2 



2mJ) 



(B-1) 



SS WW, ZZ, ff all proceed via s-channel Higgs exchange. 
The corresponding (ovrei) are; 



B. RG equations for and 

The 1-loop RG equations for the non-minimal gravity cou- 
plings and t,s are obtained as follows. For a general theory 
of scalars (|), with mass terms and non-minimal couplings in 
the Lagrangian 



L D ^miji^ii^j + ^^ijdfidfj , 
the 1-loop bare and renormaUzed / are related by ||43| 
ij = (^ki- -^hi j Z2ij + -5,7 , 



(A- 18) 



(A-19) 




2 \ 1/2 



%%m}.((Ami.-m}? + mlTV 



{OV,-el)zZ = 



(B-2) 




(B-3) 



where Z^^- ■ is the mass renormalization, 



- 7'" 



'2ij"'kl ■ 



(A-20) 



Therefore the RG equations for t,ij are related to the mass 
anomalous dimensions 7^' by 



Im 1 1 



(A-21) 



Y jj can be easily derived by applying the scalar potential RG 
equations to the 1-loop effectivepotential in order to obtain 
the p-function of the mass term [26], p„,2 ,y = 7^f,y'Wai- 

Finally, the RG equations must be modified at large s or 
h by suppressing the propagator for the corresponding real 
scalar field. Note that one does not suppress all the compo- 
nents of the Higgs doublet for the h direction nor the complex 
component of S in the s direction. 

Appendix B: Gauge Singlet Scalar Dark Matter Density 



In this appendix we give the S annihilation cross-section 
times relative velocity, (avrei), and the resulting dark matter 



and 



ff ' 



1 



3/2 



(4m2-m2) +mlTl) \ 

(B-4) 

Here the fermion Yukawa coupling is Xf = mf/v where v — 
246.22 GeV and mf is the fermion mass. Fh is the Higgs decay 
width. (Fermions should be summed over colours.) 

The dark matter density is calculated using the Lee- 
Weinberg approximation. For real S the present total mass 
density in S scalars is jioll 



Pc 



8{Ty) 



K 



(1-3W2) 



g{Tfs) TyXfsiOa 




(B-5) 
and 



where Tfs is the S freeze-out temperature, Xfs ^ Tfs/nis 

K = {ATi^g{Tfs)/A5Ml)^'^ where g{T) is the effective num- 
ber of relativistic degrees of freedom. The density for com- 
plex S is twice that for real S, due to the additional degree of 
freedom. 
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